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ABSTRACT 

Whithin the framework of the Planck satellite polarisation calibration, we present a study 
of the Crab Nebula spectral energy distribution (SED) over more than 6 decades in frequency 
ranging from 1 to 10 6 GHz (from 299 to 2.99 x 10~ 4 mm). The Planck satellite mission observes 
the sky from 30 to 857 GHz (from 9.99 to 0.3498 mm) and therefore we focus on the millimetre 
region. We use radio and submillimetre data from the WMAP satellite between 23 and 94 GHz 
(from 13 to 3.18 mm) and from the Archeops balloon experiment between 143 (2.1 mm) and 
545 GHz (0.55 mm), and a compendium of other Crab Nebula observations. The Crab SED is 
compared to models including three main components : synchrotron which is responsible for the 
emission at low and at high frequencies, dust which explains the excess of flux observed by the 
IRAS satellite and an extra component on the millimetre regime. From this analysis we conclude 
that the unpolarised emission of the Crab Nebula at microwave and millimetre wavelengths is the 
same synchrotron emission that the one observed in the radio domain. Therefore, we expect the 
millimetre emission of the Crab nebula to be polarised with the same degree of polarisation and 
orientation than the radio emission. We set upper limits on the possible errors induced by any 
millimetre extra component on the reconstruction of the degree and angle of polarisation at the 
percent level as a maximum. This result strongly supports the choice by the Planck collaboration 
of the Crab nebula emission for performing polarisation cross-checks in the range 30 (299 mm) 
to 353 GHz (0.849 mm). 

Subject headings: Cosmic microwave background - Cosmology: observations, calibration 



1. Introduction 

As the strongest source of synchrotron radia- 
tion in our Galaxy, the pulsar-powered Crab neb- 
ula ( Taurus A) is a well studied astrophysical 
object and it is therefore used for calibration pur- 
pose. This will be the case for the Planck satellite 
mission which will use the Crab Nebula for polar- 
isation cross-checks in the frequency range from 
30 to 353 GHz. A good understanding of the SED 
of the source as well as of the total intensity flux 
within the Planck beam will be required for an ac- 
curate determination of the angle of polarisation 



of the detectors and of a possible cross polarisation 
effect between detector as they limit the accuracy 
to which the CMB polarised angular power spec- 
tra will be measured. 

The emission spectrum of the A.D. 1054 su- 
pernova remnant has been the subject of a host 
of investigations over several decades in fre- 
quency. The radio spectrum is known to exhibit 
a synchrotron power l aw with a spectr al index 
/3 ~ -0.299 ± 0.009 (|Baars et al 1977T ). This 
continuum from radio synchrotron seems to be 
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fading with a rate a = (-0.167 ± 0.015)%yr" 1 
( Aller fc Reynolds 19851 ). At higher frequency, 
above 10 4 GHz , the observation are also consis- 
tent with synchrotron emission with a power-law 
of spectral index -0.73 dVeron-Cettv fc Woltier 1993h . 
The data from IRAS sate llite dMarsden et al. 1984 ) 
have been reanalyzed by (jStrom fc Greidanus 1992 ) 
revealing a significant excess of emission over the 
low frequency synchrotron spectrum, well ex- 
plained by a single dust component at a tem- 
perature T ~ 46 K, thus requiring a 0.02 M 
dust mass. Using MPifrR bolometer arrays at 
the IRAM 30 m telescope, (|Bandiera et al. 20021 ) 
gave the first evidence for a new component at 
milllimetre wavelengths. They have shown that 
this 1.3 mm (230 GHz) excess flux cannot be 
interpreted as emission from a dust component 
whereas the data may be consistent with a low 
energy cutoff in the energy distribution of the 
emitting particles. However, at 847/xm (353 GHz) 
(jGreen et al. 20041 ) found a good agreement with 
the canonical radio synchrotron emission and 
hence no need for an extra component. As the 
Planck satellite mission will use the Crab nebula 
emission to perform polarisation cross-check in 
the range 30 to 353 GHz, a deep knowledge of the 
physical origin of this emission is needed. 



For this purpose, we use observations of the 
Crab Nebula by the W MAP satellite a t 23, 
33, 41, 61 and 94 GHz (|Page et al. 20071 ) and 
by the Archeops ball oon experiment at 143 , 
217, 345 and 545 GHz (jMacias-Perez et al. 20071 : 
iDesert et al. 20081 ) . These data in addition to al- 
ready published data are used to study the Crab 
Nebula SED. The paper is organized as follows. In 
Sect. [5] we present the re-analysis of the Archeops 
and WMAP data. Sect. [3] presents the SED of 
the Crab Nebula from 1 to 10 6 GHz and compares 
it to a model including the synchrotron and dust 
well known components. In Sect. 0] we perform a 
coherent analysis of the Crab SED over the full fre- 
quency range adding an extra component to the 
previous model to account for the possible mil- 
limetre excess. The implications of our results for 
the Planck polarisation calibration are discussed 
in Sect. [5J Summary and conclusions are given in 
Sect.© 



Re-analysis of the Archeops and WMAP 
observations 




Fig. 1. — Integrated flux at 23 GHz as a func- 
tion of the angular distance to the Crab nebula 
center before (black) and after (red) background 
subtraction. 



In this paper we are mainly interested in the 
microwave and millimetre regime for which the 
best currently a vailable data are the WMAP 
(jPaee et al. 2007h and Archeops (|Desert et al. 20081 ) 
data. To avoid any bias on our analysis by a 
difference of treatment of these two data sets 
we have decided to re-estimate the Crab nebula 
unpolarised emission from 23 t o 545 GHz using 

the publicly availab le WMAP (IPaee et al. 20071: 

iHinshaw et al. 20091 ) and Archeops ([Marias-Perez et al. 20071 ) 
intensity maps. Furthermore, to make the full 
data set as homogeneous as possible, the total 
intensity flux for each frequency was estimated 
using standard radial photometry techniques to 
mimic the flux integration performed for the high 
resolution data sets. Notice that we integrate far 
from the radio extension of the Crab nebula, ~ 
5 arcminutes radius, to recover the Crab nebula 
flux diluted over the entire beam pattern, There- 
fore, a key issue for this analysis is the subtrac- 
tion of the background galactic emission which is 
mainly due to synchrotron and free-free emissions 
at low frequencies and dust at high frequencies. 
As the WMAP and Archeops maps are available 
on Healpix format we have develop our own radial 
photometry software for this pixelization scheme. 



The results of this analysis are presented on ta- 
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ble [TJ Fo r Archeops we obta in similar values than 
those in ( Desert et al. 20081 ) with slightly different 
error bars. However, for WMAP we obtain values 
8 % larger than those in (|Page et al. 20071 ). Most 
probably the difference comes from the way the 
background subtraction is performed and the fact 
the WMAP team computed the intensity and po- 
larisation simultaneously. It is important to stress 
that the difference with respect to the WMAP 
team is within the 1-er error bars at high frequency 
and goes up to 3-cr at 23 GHz. To illustrate the 
problem of background subraction we show on fig- 
ure[T]the radial profile centered at the Crab nebula 
position before (black) and after (red) background 
subtraction at 23 GHz. For this paper we use our 
own estimation of the Crab nebula flux in inten- 
sity but we have noticed that the final conclusions 
are not changed significantly by the choice of data 
set. 




Frequency (GHz) 



Fig. 2.— SED of the Crab Nebula from 1 to 
100 GHz. Data samples are represented together, 
with he best-fit power law model to the data (see 
text for details). 



3. SED of the Crab Nebula 

We present in this section a coherent analysis of 
the Crab SED in the range 1 to 10 6 GHz based on 
a compendium of observations shown in Table [1] 
Notice that to be able to directly compare to the 
Archeops and WMAP data, we have chosen only 
those data sets for which integrated fluxes over the 
full extension of the Crab Nebula are available. 

3.1. Low- frequency synchrotron emission 

An accurate determination of the low-frequency 
synchrotron component is necessary to assess the 
synchrotron contribution at mm frequencies. Any 
inter-comparison of low-frequency radio obser- 
vations of the Crab nebula must take into ac- 
count it s well-known secular d ecrease. In par- 
ticular (|Aller fc Reynolds 19851 ) have estimated 
a secular decrease in the flux density at a rate 
a — (—0.167 ± 0.015) %yr~ x from observations at 
8 GHz over the period 1968 to 1984. This result 
is in good agreement with other studies at lower 
frequencies: for example a = (—0.18 ±0.1) %yr~ Y 
over the period 1977 to 2000 at 927 MHz by 
(IVinvaikin 2005h . 

All these measurements are in fair agree- 
ment with theoretical evaluations of the evo- 
l ution of pulsar driven sup ernova remnants by 
( Reynolds fc Chevalier 19841 ) which predicts a 



For this paper the value a = —0.167 %yr~ x is 
chosen for the fading of the Crab Nebula and 
all data are converted to a common observation 
date, 01/01/2003. In figure H] we trace the flux 
of the Crab Nebula as a function of frequency for 
the fading corrected low- frequency data in table [T] 
ranging from 1 to 143 GHz. 

We observe a large decrease of flux with increas- 
ing frequency which can be represented by a power 
law of the form Ai( 1( j / Hz ) /31 . By \ 2 minimiza- 
tion, we obtain for the low-frequency data up to 
100 GHz 
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-0.296 ± 0.006 and A 1 = (973 ± 19) Jy 



with x 2 /N(j f = 1.03 and this model is shown as 
a solid line on Fig [5] This result is in agreement 
with the "canonical" v alue f3 ~ -0.299 ± 0.009 
from dBaars et al. 19771 ) at the l-a level. 



3.2. High frequency synchrotron and dust 



emission 



The Crab synchrotron emission described above 
evolves at higher frequencies, around 10 4 GHz, to- 
wards a much harder SED with a spectral index of 

0.73 (see for example (|Veron-Cettv fc Woltier 19931 ) ). 

To accurately estimate the synchrotron emission 



ranging from —0.16% to — 0.4%yr 



1 l lKovalenko et aL 19941) obtained /3i = (-0.27 ± 0.04) 



3 



v (GHz) 


S„ (Jy) 


AS,, (Jy) 


Central Epoch 


Reference 


1.117 


990.0 


59.4 


1969.9 


fVinoaradova et al. 1971") 


1.304 


980.0 


58.8 


1969.9 


fVinosradova et al. 1971) 


1.4 


930.0 


46.5 


1963 


fKellcrmann et al. 1969) 


1.765 


940.0 


56.4 


1969.9 


fVinoaradova et al. 1971) 


2.0 


840.0 


50.4 


1969.3 


fDmitrenko et al. 1970) 


2.29 


810.0 


48.6 


1969.3 


fDmitrenko et al. 1970) 


2.74 


795.0 


47.7 


1969.3 


fDmitrenko et al. 1970) 


3.15 


700.0 


24.5 


1964.4 


fMedd 1972) 


3.38 


718.0 


43.1 


1969.3 


fDmitrenko et al. 1970) 


3.96 


646.0 


38.8 


1969.3 


fDmitrenko et al. 1970) 


4.08 


687.0 


20.6 


1964.8 


fPenzias & Wilson 1965) 


5.0 


680 


34 


1963 


fKellcrmann et al. 1969) 


6.66 


577.2 


20.2 


1965. 


fMedd 1972) 


8.25 


563.0 


22.5 


1965.9 


fAllen & Barrett 1967) 


13.49 


524.0 


19.9 


1969.9 


fMedd 1972) 


15.5 


461 


24 


1965.9 


fAllen & Barrett 1967) 


16.0 


447.0 


15.6 


1970.6 


fWrixon et al. 1972) 


22.285 


397 


16.0 


1973.1 


f Janssen et al. 1974) 


22.5 


395 


7 


2003 


This paper 


31.41 


387 


72 


1966.7 


fHobbs et al. 1968) 


32.8 


340 


5 


2003 


This paper 


34.9 


340 


68 


1967.3 


fKalatrhan & Wulfsbere 1967) 


40.4 


323 


8.0 


2003 


This paper 


60.2 


294 


10.0 


2003 


This paper 


92.9 


285 


16.0 


2003 


This paper 


111,1 


290 


35 


1973.5 


(Zabolotnvi et al. 1976) 


143 


231 


32 


2002 


This paper 


217 


182 


38 


2002 


This paper 


230 


260 


52 


2000 


fBandiera et al. 2002) 


250 


204 


32 


1985.3 


fMezger et al. 1986) 


300 


194.0 


19.4 


1983 


fChini et al. 1984) 


300 


131 


42 


1978.75 


fWriffht et al. 1979) 


300 


300 


80 


1976.0 


fWerner et al. 1977) 


347 


190 


19 


1999.8 


f Green et al. 2004) 


353 


186 


34 


2002 


This paper 


545 


237 


68 


2002 


This paper 


750 


158 


63 


1978.75 


fWriffht et al. 1979) 


1000 


135 


41 


1978.75 


fWriffht et al. 1979) 


3000 


184 


13 


1983.5 


f Strom & Greidanus 1992) 


5000 


210 


8 


1983.5 


f Strom & Greidanus 1992) 


12 x 10 3 


67 


4 


1983.5 


f Strom & Greidanus 1992) 


25 x 10 3 


37 


1 


1983.5 


f Strom & Greidanus 1992) 


3.246 x 10 5 


6.57 


0.66 


1989 


fVeron-Cettv & Woltier 1993) 


4.651 x 10 5 


4.78 


0.48 


1989 


fVeron-Cettv & Woltier 1993) 


5.593 x 10 5 


4.23 


0.42 


1989 


fVeron-Cettv & Woltier 1993) 


7.878 x 10 5 


3.22 


0.32 


1989 


fVeron-Cettv & Woltier 1993) 



Table 1: Compendium of Crab Nebula observations from 1 to 10 6 GHz. Fluxes (S„) are presented in Jy. 
For Veron-Cetty & Woltjer 1993 and Chini et al. 1984 a conservative 10% error has been chosen to account 
for extrapolation errors. The central epoch of observation is also indicated. This is used for the evaluation 
of the fading effect of the low frequency synchrotron component, up to 100 GHz. Data values labeled This 
paper for Archeops and WMAP are revaluated using |he method described in Section 2. 




Fig. 3. — High frequency SED of the Crab Nebula. 
Data samples are represented together with the. 
best-fit power law model to the data (see text for 
details). 



properties at high frequency we have fitted the 
data from 10 4 to 10 6 GHz presented in table Q] 
to a power law of the form A 2 ^' 32 . By x 2 min- 
imization, we found the data to be well fitted, 
X 2 /Ndof = 0.155, by a power law of parameters 

/3 2 = (-0.698 ±0.018) and A 2 = (43.5±8.8) x 10 3 Jy 

Figure [3] represents the high-frequency data and 
the best-fit power law model is overplotted. 

Finally, the infrared satellite observatory IRAS 
has revealed significant excess emission above 

this synchrotron co ntribution aro und 50 ^im 

( Marsden et al. 19841 ). As shown by ( Strom fc Greidanus 



Fig. 4. — Residual to the standard model of the 
Crab nebula emission assuming a synchrotron and 
a dust component, in the range 10 to 2000 GHz 
(see text for details). Notice that the millimetric 
data were not used in the fit. We show the full 
data set in the region of interest. 



compare it to the data in table [TJ Fig 2] shows 
the residuals to the canonical model from 10 to 
2 x 10 4 GHz. We observe that there is no signifi- 
cant excess of power in the millimetric regime from 
100 to 1000 GHz except for the 545 GHz Archeops 
data sample which presents only a 1.5 sigma ex- 
cess. Indeed, the x 2 /Ndof for the null hypothesis 
is 0.69. 

4. Refined modeling 

In the previous section we have proved that the 
3ra$2i)netric data in table [H from 100 to 1000 GHz, 



this can be explained, after careful removing of the 
synchrotron component, by a single dust compo- 
nent described by a modified black body of emis- 
sivity f3 = 2 at T= 46 ±3 K, requiring a dust mass 
of 0.02 M©. 

3.3. Millimetric excess 

To evaluate a possible millimetric excess of flux 
in the range 100 to 1000 GHz we assume the above 
canonical modeling of the Crab nebula SED: a syn- 
chrotron component with a spectral index break at 
high frequency and a dust component at infrared 
wavelengths. We can thus reconstruct the Crab 
nebula emission at the millimetric frequencies and 



are compatible with the canonical model assuming 
single synchrotron and dust components. How- 
ever, it is interesting to check if an extra compo- 
nent may i mprove significantly t he fit to the data. 
Following (jBandiera et al. 20021 ) we consider ei- 
ther an extra low-temperature dust emission or 
an extra synchrotron component. Thus, the three 
component model is defined as follows 

1. Canonical synchrotron that is described by 
four parameters: spectral index and ampli- 
tude for the low and high frequency emis- 
sion. At high frequency we consider both the 
amplitude and the spectral index fixed and 
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set them to the values obtained in Sect. 13.21 
At low frequency we fix the spectral index 
to the value obtained in Sect. 13.11 and the 
amplitude, As, is fitted. We also assume a 
constant fading of a — —0.167 %yr~~ x as be- 
fore. 



2. Canonical dust (following (|Strom fc Greidanus 
described by a modified black body with 
three free parameters, Ajj, Tr> and /3d that 
represent the amplitude, temperature and 
spectral index respectively. 

3. One of the extra components as described 
below. 

4.1. Extra low-temperature dust emission 



In Sect. 13.31 we concluded that at 545 GHz the 
difference between the standard model and the 
data is at its maximum. In the case of an extra 
dust component this implies very low temperature 
dust (LTD) in the range from 4 to 10 K. To model 
this component we assume a modified black body 
spectrum with three free parameters Altd, Tltd 
and (3ltd that represent the amplitude, tempera- 
ture and spectral index respectively. The best-fit 
model to the data is found by x 2 minimization on 
the full frequency range from 1 to 10 6 GHz. 

Figure [5] presents on the left panel the observa- 
tional data (black) and the global best-fit model 
to the data (in red). The canonical synchrotron 
component is shown on light blue and the canon- 
ical dust in orange. The extra low temperature 
dust component is represented in blue. The right 
panel presents the residuals to the best fit-model 
in the frequency range of interest. The parameters 
and error bars for the best-fit model to the data 
and the percentage of flux at 545 GHz associated 
to the extra component with respect to the to- 
tal flux are given in table El We also present the 
X 2 /Ndof values for the global fit in the frequency 
range from 1 to 10 6 GHz and on the millimetric 
range from 100 to 1000 GHz. 

We obtain a good global fit to the data as 
shown by the \ 2 l^dof value. The best-fit pa- 
rameters obtained for the canonical synchrotron 
and dust components are in good agreement with 
those presented in Sect. [31 Comparing with 



(jStrom fc Greidanus 19921 1 . we found the same 
dust temperature, 46 ± 1 K, with an error bar im- 
proved by a factor of three, as we carefully account 
for the canonical synchrotron spectrum. The am- 
plitude of the extra component is compatible with 
zero and therefore we conclude that there are no 
— i evidence for an extra component in the form of 
^ fow-temperature dust. Furthermore, we observe 
that on the one hand the data require extremely 
low temperatures of 5 K with rather large error 
bars, therefore dust masses of ~ 230 M Q , making 
the model rather unrealistic. On the other hand, 
the improvement of the X 2 /Ndof in the millimetric 
region between 100 and 1000 GHz is not signif- 
icant to justify the addition of the three extra 
parameters required by the low-temperature dust 
component. This is also clear on the right panel 
of figure [5] where we represent the residuals to 
the best-fit model to the data on the millimetric 
regime. 

4.2. Extra synchrotron component 

For the e xtra synchrotron com ponent we con- 
sider, as in ([Bandiera et al. 20021 ). that the dis- 
tribution of energy of the relativistic electrons 
responsible for the emission is well represented by 
a power law with spectral index in the range 1 
to 3 and present a low-energy cutoff. To account 
for an excess of flux in the millimeter regime, the 
critical frequency corresponding to the lowest en- 
ergy electrons must be somewhere in the range 
200 to 600 GHz. In total the low-energy cutoff 
synchrotron (LECS) model has three parameters, 
p the spectral index of the electron energy distri- 
bution, v c the low-energy cutoff critical frequency 
and Alecs, & normalization coefficient. The best- 
fit to the data is found by x 2 minimization. 

The parameters and error bars of the best-fit to 
the data for the three data sets as well as the per- 
centage of flux due to the extra component with 
respect to the total flux at 545 GHz are given on 
table [3l We also present X 2 /Ndof values for the 
best-fit to the data on the full data set from 1 to 
10 6 GHz and on the millimetric regime from 100 
to 1000 GHz. The left panel of figure [5] presents 
the data in black and the global best-fit model to 
the data is represented in red. The canonical syn- 
chrotron component is shown on light blue and the 
canonical dust in orange. The extra low-energy 
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Fig. 5. — Refined modelling of the SED of the Crab nebula in the frequency range from 1 to 10 6 GHz 
assuming an extra dust component at low-temperature. Left panel: We represent the data from table Q] 
(black), the best-fit model to the data (red) and the synchrotron (light blue), the main dust (orange) and 
the extra dust (dark blue) components associated to it. Right panel: residuals in the millimetric regime for 
the top panel models. 



X 2 /Ndof 


As (Jy) 


A D (Jy) 


Td (K) 


Pd 


Altd (Jy) 


Tltd (K) 


Pltd 


LTD % 


Xmm/Ndof 


0.73 


971 ±9 


128 ±8 


45.9 ± 1.2 


1.93 ±0.50 


32 ±24 


5 ±23 


3.6 ±2.4 


26 ±20 


0.88 



Table 2: Best-fit model parameters and errors for the extra dust component model. \ 2 l^dof for the full 
data set and on the millimetric regime from 100 to 1000 GHz. LTD % stands for the percentage of flux at 
545 GHz due to the extra dust component with respect to the total flux. 



X Z /N dof 


A s (Jy) 


A D (Jy) 


Td (K) 


Pd 


Alecs (Jy) 


v c (GHz) 


P 


LECS % 


Xmm/Ndof 


0.76 


965 ± 14 


129 ±9 


46.0 ± 1.0 


1.89 ±0.4 


5±6 


470 ±317 


3 ±12 


3±3 


0.93 



Table 3: Best-fit model parameters and errors for the extra syncrotron component model. \ 2 /^dof for the 
full data set and on the millimetric regime from 100 to 1000 GHz. LECS % stands for the percentage of flux 
at 545 GHz due to the extra synchrotron component with respect to the total flux. 



7 



10° 10' 10 2 10° 10* 10 s 10' 

Frequency {GHz) 



100 1000 10000 

Frequency (GHz) 



Fig. 6. — Refined modelling of the SED of the Crab nebula in the frequency range from 1 to 10 6 GHz 
assuming an extra synchrotron component. Left panel: We represent the data from tabled] (black), the 
best-fit model to the data (red) and the main synchrotron (light blue), the extra synchrotron (dark blue) 
and the dust components (orange) associated to it. Right panel: residuals in the millimetric regime for the 
top panel models. 



cutoff synchrotron component is represented in 
blue. 

As above, we obtain a good global fit to the 
data as shown by the \ 2 l^dof value. The best- 
fit parameters obtained for the canonical syn- 
chrotron and dust components are in good agree- 
ment with those present ed in Sect. |3] and on 
(jStrom fc Greidanus 19921 ) . The amplitude of the 
extra synchrotron component is compatible with 
zero as well as the value of the spectral index of the 
electron energy distribution. We therefore con- 
clude that there is no evidence of an extra compo- 
nent in the form of low-energy cutoff synchrotron. 
In addition, we observe that the x 2 l^dof in the 
millimetric region from 100 to 1000 GHz is worse 
than the one obtained in Sect 13.31 assuming the 
canonical model only. The right panel of figure [5] 
represents the residuals to the best-fit model to 
the data on the millimetric regime. From this we 
conclude that the fit to the data in this frequency 
regime is not improved by adding an extra syn- 
chrotron component. 



We have also performed the analysis of the data 
set considering the p and v r parameters fixed an d 
set to the values quoted by ( Bandiera et al. 20021 ) . 
The analysis shows that the amplitude of the extra 
synchrotron component is compatible with zero 
and therefore we conclude that the data show no 
evidence of an extra synchrotron component. 

5. Implications for the Planck satellite po- 
larisation calibration 

From the previous section, we conclude that the 
intensity emission of the Crab nebula at millime- 
tre and submillimetre wavelengths is dominated 
by the well known synchrotron radiation observed 
at radio wavelengths. The data show no evidence 
of an extra synchrotron component and evidence 
for an extra dust component is not significant. 
Therefore, we can postulate that the Crab nebula 
emission at millimetre and submillimetre wave- 
lengths is produced by the same relativistic elec- 
trons that the ones producing the radio emission. 
Then it is reasonable to guess that Crab nebula 
presents the same polarisation properties at radio 
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100 200 300 400 50 100 150 50 100 150 

Frequency (Hz) Polarisation angle (degrees) Polarisation angle (degrees) 

Fig. 7. — Errors induced in the degree and angle of polarisation by an extra low-temperature dust component. 
Left plot: relative error on the degree of polarisation on the case of an unpolarised extra component as a 
function of the frequency of the polarised Planck channels. Middle plot and right plot: relative error on the 
degree of polarisation and polarisation angle error on the case of an extra polarised component as a function 
of the polarisation angle of this component. The dashed and solid lines correspond to an extra component 
with a degree of polarisation of 5 % and 10 % respectively. From blue to red we represent the values for the 
different polarised Planck channels from 30 to 353 GHz. 




Frequency (Hz) Polarisation angle (degrees) Polarisation angle (degrees) 

Fig. 8. — Errors induced in the degree and angle of polarisation by an extra low-energy cutoff synchrotron 
component. Left plot: relative error on the degree of polarisation on the case of an unpolarised extra 
component as a function of the frequency of the polarised Planck channels. Middle plot and right plot: 
relative error on the degree of polarisation and polarisation angle error on the case of an extra polarised 
component as a function of the polarisation angle of this component. From blue to red we represent the 
values for the different polarised Planck channels from 30 to 353 GHz. 
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and millimetre wavelengths. These conclusions are 
also supported by the Crab polarised emission ob- 
served from 23 to 94 GHz by the WMAP satellite 
( Page et al. 20071 ) when compared to the 363 GHz 
SCU BA measurements (Greaves et al. 20031 ) and 



the (|Flett fc Murray 199J ) 273 GHz data. From 
these, it makes sense to use the low-frequency 
observations of the Crab nebula to cross-check 
and eventually update the knowledge of the po- 
larisation characteristics of the high-frequency 
instruments as proposed by the Planck col- 
laboration. As a summary and using the re- 
sults from dPaee et al. 20071: iGrea ves et al. 2051 : 
Flett fc Murray 19911 ) we expect the Crab nebula 
emission at the Planck frequencies to be polarised 
with a degree of polarisation of 8 to 9 % and a po- 
larisation angle of 150° in equatorial coordinates 
(|Aumont 20091 ). Then, at 30 GHz and 353 GHz 
we expect the total intensity to be 351 ± 7 and 
173 ± 5 Jy and the polarised intensity about 28 
and 15 Jy, respectively. 

However, we can estimate upper limits on the 
errors induced on the determination of the degree 
of polarisation and polarisation angle of the Crab 
nebula emission at the Planck observation frequen- 
cies by the presence of an extra component. For 
this we use the results obtained on the previous 
section and we assume two different cases: po- 
larised and unpolarised extra component. Fig- 
ure [7] shows 1-cr upper limits on the case of a low- 
temperature dust extra component. For the left 
plot we assume an unpolarised extra component 
and we present the error on the determination of 
the degree of polarisation as a function of the ob- 
servation frequency. This error is below 1 % for 
the Planck CMB channels from 70 to 217 GHz 
and rise up to 8% at 353 GHz. The middle and 
right plot assume a polarised extra component and 
show the error on the polarisation angle and de- 
gree of polarisation, respectively, as a function of 
the polarisation angle of the extra component with 
respect to the main synchrotron component. From 
blue to red we show the errors from 30 to 353 GHz. 
The dotted-dashed and solid curve assume a de- 
gree of polarisation of 5 and 10 % for the extra 
component. As before the error on the degree of 
polarisation is below 1 % at the CMB channels 
and from 12 to 18 % at 353 GHz. For the polar- 
isation angle the error is below 0.2° for the CMB 



channels and at maximum of 3° at 353 GHz. No- 
tice that the increase in the errors at 353 GHz is 
mainly due to the large value and large error bars 
of the Archeops data at 545 GHz and we do not 
think it is significant. Figure [5] shows the errors 
induced by an extra low-energy cutoff synchrotron 
component either unpolarised or polarised with a 
degree of polarisation equivalent to the one ex- 
pected for the canonical synchrotron component. 
In this case the errors are much lower being well 
below 1 % and 0.2° for the degree and angle of 
polarisation, respectively, at all polarised Planck 
frequencies. 

6. Summary and conclusions 

Within the framework of the polarisation cali- 
bration of the Planck satellite emission we present 
in this paper a global analysis of the SED of the 
Crab Nebula in the frequency range from 1 to 
10 6 GHz. For this purpose we have used new data 
from the WMAP satellite dPaee et al 20071) and 
from the Archeops balloon experiment ( Desert et al. 2008 
in addition to data currently available. We fo- 
cus on the centimetric and millimetric regime as 
Planck observes the sky from 30 to 857 GHz. 

We have first shown that the data are com- 
patible with the canonical model of Crab neb- 
ula emission which assumes a synchrotron com- 
ponent at low and high frequencies, plus a dust 
component at the micrometre wavelengths. To 
check if an extra component may improve signif- 
icantly the fit to the data, we consider cither an 
extra low-temperature dust emission or an extra 
low-energy cutoff syn chrotron component, as in 
(jBandiera et al. 20021 ) . In both cases, we conclude 
that the current data present no evidence of an ex- 
tra component. Therefore, the un-polarised mil- 
limetre flux of the Crab nebula is well represented 
by the following synchrotron power-law model 



F v = 973 ± 19 



Jy( r 



GHz' 



(-0.296 ± 0.006) 

exp (a(T obs 
(1) 

1 is the synchrotron 



2003)) 



where a = — 1.67x 10 yr 
fading and v and T Q b s are the frequency and the 
date of observation in years A.C., respectively. 

From above, we can conclude that the millime- 
tre emission of the Crab nebula has the same phys- 
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ical origin than the radio synchrotron emission and 
therefore it is expected to be polarised with the 
same degree of polarisation and the same orien- 
tation. From the current data set we estimate 
that the errors on the reconstruction of the degree 
and angle of polarisation on the Planck cosmolog- 
ical channels induced by an extra component will 
be well below 1 %. This result strongly supports 
the choice by the Planck collaboration of the Crab 
nebula emission for performing po larisation cross - 
checks in the range 30 to 353 GHz (lAumont 2009h . 
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